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APPIJCATION  OF  THE  RATE  DIAGRAM  TECHNIQUE  TO  THE 
ANALYSIS  AND  DESIGN  OF  SPACE  VEHICLE  ON-OFF 
ATTITUDE  CONTROL  SYSTEMS 

* 

H.  Patapoff 

ABSTRACT 

Preliminary  single  axis  analysis  and  design  of  a  space  vehicle  control  system, 
consisting  of  a  torque  producing  device  and  a  dead -bard  id. thin  lAiich  there 
are  no  exteraally  or  intomal.ly  applied  torques,  can  be  accomplished  quite 
effectively  by  a  simple  graphical  technique  termed  the  "Rate  Diagram"  method. 
Two  characteristics  of  typica]  space  vehicle  attitude  control  problems  vAich 
make  this  method  effective  are  the  undamped  rigid  body  motion  of  the  vehicle 
and  the  associated  low  angular  rfites.  These  charucteristics  satisfy  the  two 
basic  assumptions  required  for  the  use  of  this  technique;  namely,  that  the 
vehicle  angular  acceleration  is  proportional  to  the  applied  torque,  and  that 
system  transients  effectively  decay  prior  to  control  torque  application.  The 
system,  exclusive  of  the  vehicle  dynamics,  can  be  of  any  order,  thus  allow¬ 
ing  sensor  dynamics,  shaping  networks,  control  torque  characteristics,  etc., 
to  be  included. 

A  "Rate  Diagram"  is  5lmp.l.y  a  plot  of  the  vehicle  angular  rate  at  control 
torque  removal  versus  tlie  rate  at  torcgie  application.  InJ'ormation  regarding 
system  stability,  transient  response,  and.  limit  cycle  behavior  can  be  ob- 
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tained  directly  from  such  a  diagram.  'ITiis  diagram  can  be  constructed  with 
relative  ease  through  the  use  of  the  analytical  relationships  developed  and 
presented  in  this  paper.  In  addition,  these  relationships  allow  the  analy¬ 
tical  determination  of  limit  cycle  amplitude  and  period. 

This  graphical  method  of  analysis  can  be  quite  effective  as  a  design  tool 
in  designing  systems  to  meet  specified  performance  requirements.  As  system 
parameters  or  control  torque  characteristics  are  varied,  the  resulting 
changes  in  system  behavior  can  be  noted  visually  by  the  corresponding 
changes  in  the  Rate  Diagram.  The  interj^retation  of  a  l^te  Diagram  for 
analysis  and  design  purposes  is  discussed,  and  a  detailed  treatment  of  an 
example  space  vehicle  control  system  probiem  to  exemplify  this  technique  is 
included. 


IHTRODUenON 

On-off  reaction  control  systems  have  found  wide  application  in  the  field 
of  Ejace  vehicle  attitude  control.  This  control  method  is  particularly 
attractive  because  of  its  Iniierent  simplicity  and  relatively  high  relia¬ 
bility,  being  most  effective  and  economical  in  cases  wiiei’e  only  relatively 
sma].l  correcting  torques  are  necessary  to  maintain  attitude  control  - 
especially  in  cases  where  the  allowable  attitude  excursions  are  relatively 
large . 

The  control  torque  generators  may  be  either  torque  producing  devices  which 
change  the  total  vehicle  angular  momentum  in  inertial  space,  or  momentum 
storage  devices  wlilch  transfer  angniar  momentum  between  the  vehicle  and 
rotating  Inertias  within  the  vehicle.  The  most  eoinmonly  used  torque  pro¬ 
ducing  flevlces  are  of  the  mass  ejection  type.  Compressed  gases,  llquiil 
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vaporization;  or  the  more  sophisticated,  ion-generation  are  exon^les.  The 
most  common  momentum  storage  device  is  the  veil  knovm  reaction  v^eel. 

One  of  the  basic  difficulties  associated  with  ouch  nonlinear  control  systems 
is  the  lack  of  available  analysis  or  synthesis  techniques  vdilch  are  suitable 
for  the  practicing  control  engineer.  One  approach  is  to  treat  the  problem 
as  a  conventlohal  one  applicable;  say;  to  phase  p^^ne  analysis;  ignoring 
effects  vhlch  moke  the  problem  ccjiq>llcated;  such  as  actuator  and  sensor 
dynamics.  If  there  are  no  serious  problems  associated  with  such  an  over¬ 
simplified  case;  it  la  then  hoped  that  the  results  vrill’ serve  as  a  sort 
of  first  approximation  to  the  "real"  solution. 

A  technique  is  presented  in  this  paper  in  vrfiich  the  performance  of  space 
vehicle  control  systems  can  be  attalyzed  vdlthout  the  need  of  over-simplifying 
the  control  prob].em.  Actuator  and  sensor  dynamics  and  shaping  filters 
can  be  Included  in  ouch  an  analysis.  This  technique  iS;  however;  based 
upon  two  major  assumptions;  namely,  that  a].],  ti'anoients  in  the  system 
effectively  decay  prior  to  control  torque  command,  und.  that  the  vehicle 
dynamics  can  be  represented  as  a  donb.lc  integration  of  the  control  torque. 
These  two  assumptions  have,  in  the  past,  been  valid  in  numerous  space 
vehicle  control  problems. 


GYfyrpji  DEfJCi^iprioN 

A  block  diagram  of  the  on-off  control  sysbem  considered  appears  in  Figure  1. 
I  is  the  moment  of  inertia  of  the  vehicle,  v/}iogc  angular  acceleration  is 
proportional  to  tlic  control  toi’qnc,  T^.  'J'Jic  voliiclc  attitude,  0,  is  sensed 
and  compared  vdth  the  commanded  attitude,  0^.  The  resulting  error  signal 
is  filtered  and  fed  to  the  switch.  The  switch  has  a  dead -space  of  +  0j^ 
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and  a  per-unit  hysteresis  h.  The  system  operation  can  be  explained  with 
the  aid  of  the  phase-plane  trajectory  of  Figure  2.  Assume  that  the  vehicle 
rotates  at  an  angular  rate,  and  that  the  output  signal  of  the  filter 
lies  vrlthin  the  dead  space  of  the  switch.  As  the  vehicle  rotates,  the  an¬ 
gular  error  Increases  until  the  filter  output  becomes  G^,  at  which  time 
switch  closure  occurs.  This  corresponds  to  point  1  of  Figure  2.  seconds 
after  switch  closure,  control  torque  is  applied  vdiich  tends  to  reverse  the 
direction  of  travel,  corresjxjndlng  to  point  2.  Tyie  transportation  delay 
time,  is  assumed  to  be  constant.  Control  torque  is  applied  to  the 

vehicle  until  the  output  of  the  filter  reduces  to  the  value  Gjj(l-h),  at 
vAiich  time  the  switch  opens.  This  corresponds  to  ixjint  3*  Control  torque 
is  not  removed  until  seconds  after  switch  opening,  corresponding  to 
point  4.  The  filter  output  is  again  within  the  dead-space  of  the  switch, 
and  the  vehicle  rotation  has  been  reversed  at  a  rate,  G^.  This  sequence 
is  now  repeated.  If  the  vehicle  rotation  has  not  been  reversed  after 
torque  removal,  then  torque  will  again  be  applied  until  reversal  of  rotation 
does  occur.  An  example  i^ase-plane  trajectory  for  this  occurance  is  sho  .11 
in  Figure  3. 

The  filter  is  assumed  to  have  a  rational  algebraic  transfer  function  (with 
poles  in  the  left  half  s  plane)  with  unity  zero  frequency  gain.  The  actuator 
and  sensor  dynaiDj.es  are  included  in  the  filter  transfer  function. 

The  control  torque  build-up  imd  tall-off  characteristics  are  assumed  to  be 
the  same  for  either  positive  or  negative  vehicle  accelerations.  Figure  4 
shows  typical  torque  characteristics  for  a  constcuit  torque  generator.  For 
convenience  the  control  torque  tail-off  characteristics  are  included  in  the 
time  delay,  x^.  For  the  typical  case  of  Figure  4  the  tail-off  character- 
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IstlcB  are  independent  of  the  switch-off  time  nnce  full  torque  has  been 
achieved.  If  shut-off  occurs  prior  to  full  torque  build-up, or,  in  the 
general  case,  if  tail-off  characteristics  are  a  function  of  the  control 
torque  "on-time",  the  necessary  modification  of  vs.  control  torque 
"on-time"  can  be  determined  experimentally. 

THE  RATE  DIAGRAM 


Definition 

A  rate  diagram  is  simply  a  plot  of  the  vehicle  angular  rate  at  control 

•  • 
torque  removal,  0^,  versus  the  rate  at  control  torque  application,  6^. 

Construction 

Analytical  expressions  for  use  in  constructing  such  a  diagram  are  derived 
in  the  Appendix. 

Since  the  control  system  characterlsti..s,  except  for  sign,  are  assumed 
identical  for  both  positive  and  negative  vehicle  rotations,  it  is  suff¬ 
icient  to  construct  rate  diagrams  only  for  positive  values  of  6^.  It 
is  sufficient  to  know  only  the  angular  rate  at  torque  removal,  6^,  re- 

lative  to  0  . 

-  0 

Due  to  the  fact  that  the  control  torque  tends  to  reverse  the  direction  of 
travel,  the  vehicle  final  angular  rate,  6^,  for  an  infinitesimally  small 
positive  initial  rate,  0^,  is  negative.  Also,  if  the  final  rate  has  the 
same  sign  as  the  initial  rate,  its  magnitude  will  be  less  than  that  of 

the  initial  rate.  This,  in  effect,  means  that  the  rate  curve  (for  posl- 

.  •  •  * 

tlve  0^)  always  lies  below  the  line  0^  -  0^  =  0. 
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Inte  rpre tat Ion 

1.  Stability 

The  stability  of  the  control  system  can  be  determined  by  the  addition 
of  the  line,  (which  will  be  termed  the  "minus-one"  line) 

to  the  rate  diagram.  If  the  rate  curve  lies  eveiyvdiere  above  this  line, 
the  rate  amplitude  continually  decreases;  indicating  a  stable  system. 

If  the  rate  curve  lies  everyvdiere  below  this  line,  then  the  rate 
amplitude  continually  increases,  indicating  an  unstable  system. 

Example  diagrams  appear  in  Figure 

In  general,  the  "minus-one"  line  (©^  +  ©^  =  O)  will  Intersect  the 
rate  curve,  indicating  regions  of  stability  and  instability.  All 
points  of  intersection  indicate  the  existence  of  syranetrlcal  limit 
cycles,  since  ©^  =  -  0^*  Whether  these  limit  cycles  arc  stable  or 
unstable  depends  upon  the  manner  in  vrtiich  the  rate  curve  inter¬ 
sects  the  minus-one  line.  Unsymroetrical  limit  cycles  (\rtiere  the 
magnitudes  of  the  initial  euid  final  rates  are  unequal)  can  be  found, 
if  they  exist,  by  constructing  lines  perpendicular  to  the  minus-one 
line.  If  any  of  these  lines  intersects  the  rate  curve  at  more  than 
one  point,  and  if  two  points  of  intersection  are  equidistant  from 
the  minus -one  line,  then  an  unsymmetrical  limit  cycle  is  indicated. 
Reference  to  Figure  6b  shows  that  points  1  and  2  correspond  to  an 
unsymioetrical  limit  cycle.  The  existence  and  stability  of  such  limit 
cycle  behavior  will  not  be  discussed  in  this  paper. 

2.  Synmetrical  Limit  Cycles 

Consider  the  example  rate  diagram  of  Figure  6a.  The  two  points  of 
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Intereectlon  Indicate  the  existence  of  synmetrlcal  limit  cycles. 

Points  1  and  2  represent  a  stable  and  unstable  limit  cycle,  res¬ 
pectively.  Considering  point  1,  a  sll^t  decrease  In  rate  amplitude 
results  In  a  bulld-up  of  amplitude,  and  a  slled^t  Increase  results 
In  a  decrease  In  amplitude.  The  system,  then,  undergoes  stable 
oscillations  In  the  vicinity  of  point  1.  For  point  2,  a  slight 
Increase  In  rate  amplitude  results  In  a  continual  rate  bulld-up. 
Indicating  an  unstable  system  for  rates  above  that  corresponding 
to  i>olnt  2.  A  slight  decrease  In  rate  amplitude  results  In  a  rate 
decrease  until  point  1  Is  reached. 

Ibe  cunplltudes  of  the  limit  cycle  rates  can  be  read  directly  from 
the  rate  diagram.  If ,  however,  the  limit  cycle  rate  amplitude  Is 
of  sole  Interest  (assuming  knowledge  of  existence  and  stability). 

It  can  be  detenulned  analytically  by  use  of  equations  (A-I3),  (A-15), 
and  (A-i6)  derived  In  the  Appendix. 

For  example,  assuming  a  step  function  for  the  control  torque  character- 
Istlcs,  the  limit  cycle  rate  amplitude,  0,™,  is  determined  using 
equation  (A-I5): 


-S 


\ 


®D  *  ®Lc"a  *  »LC 


28. 


where 


e  '  -  Ao 


lili 

s3 
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(2) 


7 


'8982-0005 -RU-000 


Expressing  P(8)  in  the  form: 


F(s)  = 


•  •  m-l 

A  s  +  A  ,8  + 

m  m-l 


_  n  _  n-1 

B  s  +  B  ,8  + 

n  n-1 


+  AgB  +  Aj^s  +  1 


+  BgS  '  +  Bj^s  +  1 


(3) 


or. 


F(8) 


m 

TT<^  * 

j  * 


s) 


TT  *  'p/' 

J=i 


(M 


and  expanding  into  partial  fractions  (assuming  that  the  order 

of  the  denominator  of  — is  greater  than  that  of  the  numerator): 


F(s)  1  ^2  ^1 

S  S  S  8 


J=1 


1  +  1  8 


(5) 


gives: 


-1 


F(s) 

„3 


X*  i: 

j=i 


« 


(6) 


Hie  partial  fraction  coefficients  of  equation  (5)  can  be  determined 
from  the  folloifing  relationships; 


°2  =  \  -  ®1 

Cl  =  Ag  -  Bg  -  Bi  (Aj  -  B3) 


(7) 

(8) 


8 


a5aE*.OOD5»BUiOOO 


where: 
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*1  ■  I 
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m  m 
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If  F(b)  contains  simple  poles: 
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Coniblning  equations  (l)^  {Z),  and  (6): 
/ 


LC 


T.+T„\l\ 


"I  2 


*  ^o'^R°2  "  2 


Jl,  i  K,  ^LC  V 

J=i  ^ 


\diere :  2q 

+  ^  T 

■  TT  ■ 


(17) 


Assuming  simple  poles  for  F(s),  and  substituting  equation  (l^)  Into 

(16): 


\  To2 

^o  R 


'LC 


■ 

(T.+T_ 

“ 

A  R 

°2" 

2 

X  )  (1-e  T  ) 

''o  Z_,  '  Pj  ' 


(18) 


J=1 


If  F(e)  Is  simply  a  pure  lead, 


F(s)  =  1  + 


(19) 


then  Equation  (18)  reduces  to  the  algebraic  relation; 
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4  (2  +  a^m^) 
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If  F(8)  is  a  lead-lag  filter, 
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1  +  T^S 

F(b)  =  - — 

1  +  TgS 

then  Equation  (l8)  becoineE: 


(2R) 


LC 


o  R 

Tr~ 
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+  'a”i> 


"iXo'z'"!-  'z'> 


(!-«■“) 


(23) 


where 


(24) 


a 


(25) 


Since  the  term,  <x,  in  Equation  (23)  is  a  funot.lon  of  0^^,  a  graphical 
or  iterative  solution  for  is  necessary. 

The  maximum  excursion  of  the  vehicle,  Oj^,  (see  Figure  7)  is  given  by 
the  relation; 


=:  0 
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^LC  ^LC 

-  m 


(26) 


where 


(27) 
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Here,  -ra  is  the  slope  of  the  control  torque  "on"  line.  Because  of 
the  assumption  thnt  transients  effectively  decay  prior  to  swltch-on, 
these  ’  on"  lines,  vAiich  intersect  the  6  axis  at  +  9^^,  can  be  con¬ 
structed  in  the  phase-plane.  'Fhe  construction  of  "off"  curves, 
however,  become  quite  tedious. 

Tile  period,  Tj^^,  of  the  limit  cycle  is  determined  by  noting  (see 
Figure  ?)  that  the  time  interval  between  points  1  and  2  is  simply 

®o/®LC^ 

®o  =  ®D-  -¥  (28) 


The  total  "off-time ',  per  cycle  is  then; 


T  =  li 
off  ^ 


0. 


a 


D  _  1 

m 
LC 


(29) 


wxicre 


The  time  interval  between  points  ?.  and  3  is  equal  to  t  +  t  , 

LC  R' 

^LC  switch  opening,  and  is  given  by  equation  (17)  for 

the  case  of  a  step  torque.  For  the  more  general  case  t  and  Q  are 

LC  LC 

related  by: 


(30) 


This  expression  is  obtained  from  Equation  (A-I6)  by  setting  0  =  -  ©  . 

f  o’ 

For  a  step  torque,  the  total  control  torque  "on-time"  per  cycle  is: 
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on 


‘'“lc 

JT 


(31) 


This  gives  a  limit  cycle  period  of 


LC 


’^on  ♦  Vf  ■  '* 


'D 


IJC 


6. 


-L2 

K 


1 

m 


(32) 


The  total  impulse  expended  by  the  control  nystem  per  cycle  is  equal 
to  *^10, n*  Th®  impulse  per  unit  time,  xdiich  gives  a  measure  of  fuel 

Lrt 

consumption,  is  then: 


Impulse/unit  time 


hlQ. 


LC 


(33) 


LC 


If  minimization  of  fuel  consumption  during  limit  cycling  is  the  goal 
for  a  specific  control  problem,  this  is  accomplished  analytically 
by  partial  differentiation  of  Equation  (33)  with  respect  to  the 
parameter  of  interest,  or  graphically  by  simply  plotting  Equation 
(33). 


3.  Construction  of  Fhase-Plane  Trajectories  From  the  Rate  Diagram 

In  order  to  have  a  better  understanding  of  the  interpretation  of  a 
rate  diagram,  let  us  construct  a  typical  phase-plane  trajectory 
from  the  example  rate  diagram  of  Figure  8.  Let  us  assume  a  step 
control  torque  and  an  initial  vehicle  rotation  at  an  angular  rate 
corresponding  to  point  A  of  Figure  9-  The  response  of  the  system 
for  this  Initial  rate  will  now  be  followed.  The  vehicle  rate  at 
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switch -off  (point  B)  has  n  value  that  is  opposite  in  sign  to  that 
corresponding  to  point  A.  This  means  that  the  vehicle  motion  has 
reversed  direction.  Tlie  rate  corresponding  to  B  will  now  remain 
constant  until  again  sv/ltch-on  occurs.  B,  \diich  is  no./  considered 
the  initial  rate,  gives  a  final  rate,  C,  at  switch-off.  C  is  pos¬ 
itive,  indicating  that  the  vehicle  is  still  rotating  In  the  same 
direction,  but  with  a  reduced  magnitude.  After  swltch-on  again 
occurs,  the  final  rate  becomes  that  corresponding  to  point  D  at 
switch -off.  The  vehicle  has  now  reversed  direction.  This  is  re¬ 
peated  for  point  E.  It  is  again  noted  that  control  torque  "on" 
lines  can  be  constructed  in  the  phase-plane  diagram.  These  lines 
intersect  the  0  axis  at  +  9^^,  and  Imve  a  sloixj  of  -m,  vdiere  m  is 
determined  from  equation  (27).  If  P(s)  is  expressed  in  the  form  of 
Equation  (3),  or  (A-l)  of  the  Appendix,  then  m  is  given  by: 


m 


1 

\  -  ’a 


(Sk) 


Tlie  slope  of  the  switch-on  lino  is  simp]  / 

''l  ■ 

For  this  special  caEe(step  control  torque)  control  torque  "off  curves 
can  reac’ily  be  constructed  in  the  phase-plane,  since  the  trajectories 
are  paraboiic.  However,  in  the  more  general  case,  where  the  control 
torque  is  a  function  of  time,  the  construction  of  these  "off  curves 
becomes  rather  tc(ilous  since  the  vehici.c  attitude  must  be  determined 
at  control  torque  removal .  Numerical  integration  may  be  necessary. 
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4.  Transient  Response 

Some  general  remarka  can  be  made  concerning  the  transient  behavior 
of  the  system.  If  the  rate  curve  lies  norr,  and  is  parallel,  or 
nearly  parallel,  to  the  minus-one  line,  then  convergence  is  relatively 
slow.  Rapidity  of  convergence  increases  as  the  rate  curve  nears  the 
0^  axis.  Examples  appear  in  Figure  10. 

Convergence  to  the  limit  cycle  for  rates  in  the  vicinity  of  a  stable 
limit  cycle  depends  upon  the  slope  and  sliape  of  the  rate  curve  at 
the  intersection  vdth  the  minus-one  line.  If  the  slope  of  the  rate 
curve  approaches  +  1  (perpendicular  to  the  minus-one  line)  at  the 
intersection,  convergence  becomes  clow.  Also,  if  the  slope  of  the 
curve  approaches  -  1  (parallel  to  the  minus-one  line),  convergence 
becomes  slow.  Convergence  to  the  limit  cycle  improves  as  the  slope 
of  the  rate  cuin^o  approaches  zero  at  the  point  of  intersection. 

Example  rate  diagrams  with  corresponding  example  phase -plane  tra¬ 
jectories  appear  in  Figures  11,  12,  and  13. 

DE.SIGN  APPlJCABIhlTY 

Several  basic  problems  arise  if  a  systems  approach  is  taken  in  the  design 
of  a  satellite  or  space  vehicle.  One  is  a  "proper"  choice  of  the  various 
actuating  and  sensing  combinations  wiiich  vdll  satisfy  the  performance  re¬ 
quirements  of  the  mission.  In  general  this  choice  will  not  be  unique.  In 
addition,  the  determination  of  "desired"  characteristics  of  the  various 
control  system  elements,  and  of  "proper"  si^pial  processing  confronts  the 
control  engineer.  Tlie  choice  of  a  design  criterion,  from  v/hich  various 
control  system  decifjn  configurations  can  be  evaluateii,  plagues  the  control 
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engineer  because  of  the  difficulty  of  divorcing  the  general  problem  from 
the  details  of  the  mission  and  the  operational  concept,  and  the  difficulty 
in  establishing  a  set  of  system  "values"  or  absolute  standards  which  allow 
system  "optimization".  Whether  this  set  of  values  or  standards  is  the 
"right"  one  is  a  rather  philosophical  problem  in  Itself.  Because  of  these 
basic  difficulties,  control  systems  are  not  presently  designed  using  such 
a  generalized  "rational"  synthesis  approach.  The  control  engineer  is 
forced  to  select  tentative  systems  and  conduct  a  performance  analysis  for 
each  system.  A  trade-off  study  is  then  conducted  to  compare  the  relative 
merits  of  the  tentative  systems,  and  a  final  design  is  then  chosen.  This 
final  choice  is  in  essence  an  intuitive  one.  In  which  the  designer  utilizes 
his  engineering  "Judgment'  or  ’feel". 

The  synthesis  ai>pllcabillty  of  the  rate  diagram  technique  presented  in 
this  pai)er  is  based  upon  this  design  approach.  If,  for  specific  control 
system  configurations,  the  effects  of  important  parameter  variations  upon 
system  performance  can  be  understood  with  the  aid  of  such  diagrams,  the 
control  engineer  can  at  least  develop  a  "feel  ,  however  crude  it  may  be, 
for  the  overall  design  problem.  Experience  and  engineering  Judgment  are 
still  necessary  in  the  final  system  design. 

EXAMPLE  CONTROL  SYSTEM 

TSie  effects  upon  system  performance  of  an  example  control  system  are  shown 
in  rate  diagram  form  (Figures  l4,  I6,  I8,  and  20)  for  fovu:  cases  of  para¬ 
meter  variation. 

The  ncmilnal  system  parameters  chosen  are: 
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X9 

^  B  1.0  deg/sec  (step  function  assunsd) 

B  Tp  ■  0.05  sec 

hOp  >  0*1  deg 

t.s  +  1 

K.)  -  - - — 5 - = - 

<  3-2  *  t  * 

Case  1.  -  Variation  of  Lead  Tine  Constant  (t^) 

Hie  effect  of  filter  lead  time  constant  variation  upon  the  rate  dlagran 
and  limit  cycle  an^lltude  are  shown  In  Figures  l4  emd  15 .  Hie  filter 
transfer  function  for  this  case  is: 

T-S  +  1 

F(8)  .  - 

0.1  Tj^S  +  1 

Case  2.  -  Variation  of  Time  Delay  =  Tj^) 

Hie  effect  of  time  delay  variation  upon  the  rate  diagram  and  limit  cycle 
amplitude  ere  shown  In  Figures  16  and  17 .  The  filter  transfer  function 
for  this  case  Is: 

F(.)  . 

.58  +  1 

Cwe  3.  -  Variation  of  ^steresis 

Hie  effect  of  hysteresis  variation  is  shown  in  Figures  I8  and  19 .  The 
filtet  transfer  function  is; 
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F{«)  . 

.5s  +  1 

Case  U.  -  Addition  of  Quadratic  Lag 

The  effect  of  adding  a  quadratic  lag  to  the  system  is  shown  iq  Figures 
20  and  21.  The  filter  transfer  function  is: 


F(s) 


lOs  +  1 


(s  +  1)  (  ^_  +  +  1) 

cs  2  0)  ' 

n  n 


CONCLUSIONS 

The  detailed  treatment  of  the  example  control  system  shows  the  effectiveness 
of  this  technique  in  designing  space  vehicle  on-off  control  systems.  As 
filter  parameters  or  control  system  characteristics  are  varied,  the  design 
engineer  "observes"  the  resulting  changes  in  system  behavior  by  noting  the 
corresix)nding  changes  in  the  rate  diagram.  This  can  give  a  general  "feel" 
for  system  optimization  in  meeting  the  specified  performance  requirements 
of  the  mission.  Tlie  limit  cycle  derivations  can  readily  be  applied  In  cases 
vdiere  fuel  consumption  during  limit  cycling  is  of  prime  interest. 

A  general  design  approach  is  suggested  by  specifying  the  desired  rate 
diagram  characteristics  and  deslgring  a  filter  ^diich  will  yield  the  de¬ 
sired  rate  diagram,  or  at  least  a  satisfactory  approximation.  This  filter 
synthesis  would  probably  require  computer  techniques. 
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SWITCH 


Figure  1.  System  Block  Diagram 


e 


I.  SWITCH  CLOSURE 
2  TORQUE  APPLICATION 

3.  SWITCH  OPENING 

4.  TORQUE  REMOVAL 


Figure  Z.  Example  TliPce-Plane  Trajectory 


9 


Figure  3»  Example  niaso-Plane  Trajectory 
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(a)  STABLE  (b)  UNSTABLE 


Figure  5*  li^te  Diagrams  Indicating  a  Stable  and  Unstable  System 


(o)  SYMMETRICAL  (b)  ONSYMMETRICAL 


Figure  6.  Date  Diagrams  Indicating  Limit  Cycles 
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Figure  8.  KxPinpD.c  Itote  Diagram 


Figure  9,  Typical  Phanc  Plane  Tra.lectory  for  Kxainple  System  of  Figure  8. 


(a)  SLOW 


(b)  RAPID 


Figure  10.  Example  Rate  Diagrams  Indicating  Relatively 
Slow  and  Rapid  Convergence 


Figure  11.  Example  Indicating  Relatively  Slow  Convergence 
In  Vicinity  of  Stable  Limit  Cycle 
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t 


(O)  RATE  DIAGRAM 


r 


(b)  typical  phase- plane  trajectory 


8 


Figure  12.  Example  Indicating  Relatively  Slow  Convergence 
In  Vicinity  of  Stable  Limit  Cycle 


Figure  13.  Exanqjle  Rate  Diagram  Indicating  Rapid  Convergence 
In  Vicinity  of  Stable  Limit  Cycle 
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Figure  20.  jSxample  Control  System  Rate  Diacrams  vs.  Quadratic 
lag  Natural  Frequency,  u>^  ^  S  ''2  ^ 


Figure  21.  Limit  Cycle  Hate  Amplitude,  0^^^  vp.  o,,adratic  lag  Frequency, 
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APPENDIX 

Derivation  of  Equations  for  Rate  Diagram  Construction 


Procedure 

In  constructing  a  rate  diagram  an  analytical  expression  relating  the 
vehicle  angular  rate  at  control  torque  removal  to  the  rate  at  appli¬ 
cation  is  first  derived.  This  is  accomplished  by  determining  the  "on- 
time"  of  the  control  torque,  whereby  the  change  in  vehicle  angular  rate 
can  be  determined  from  a  knowledge  of  the  torque  characteristics.  A 
plot  of  the  rate  at  torque  removal  versus  the  rate  at  torque  application 
is  the  rate  diagram. 

Assumptions 

The  following  derivations  are  based  upon  the  assumptions  that: 

1.  All  transients  in  the  system  have  effectively  decayed  prior  to 
switch -on.  In  other  words,  the  "off -time'  of  the  control  torque 
is  long  compared  to  the  system  time  constants. 

2.  The  vehicle  dynamics  can  be  represented  as  a  free  rigid  body  in 
vrtilch  the  acceleration  is  proportional  to  the  applied  torque. 

3.  The  control  system  characteristics  are  identical  for  both  positive 
or  negative  vehicle  rotations. 

For  sin^liflcatlon,  the  commanded  vehicle  attitude,  6^,  is  assumed  to  be 
zero. 
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Derivation 

Ijet  the  filter  transfer  ftinctlon,  P(s)  be  expresced  nc  the  ratio  of  the 
two  polynomials: 


F(8) 


nM  - 

D(s) 


A  a"  +  A  ,8™"^  +  .  .  .  +  A,8  +  1 
m _ m-1 _ 1 _ 

B  a”  +  B  ,b""^  +  .  .  .  +  B,b  +  1 
n  n-1  1 


(A-1) 


Since  no  torque  Is  applied  to  the  vehicle  prior  to  switch -on,  the  vehicle 
rotates  at  a  constant  rate,  6^.  The  filter  input,  then,  is  simply  a  rati^. 
The  steady-state  output  of  the  filter,  to  the  ramp  input,  6^,  Is: 


and 


ss 


0  -t 


(Ai- 


e  . 

o 


(A-2) 


(A-3) 


It  is  seen  that  the  filter  otpat  leads  the  input  by  an  amount  0  (A,-  B, ), 

O  a.  J. 

and  has  the  same  slope. 

Consider  the  case  >;here  switch-on  occurs  vdien  g  -  -  Control  torque 

application  occurs  seconds  .later,  at  v/hich  time: 


e 


(Or 


0  T.  ) 

o  A 


(A-4) 


€ 


(A-5) 
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Now  let  time  begin  at  control  torque  application.  T)ie  differential 
equation  describing  the  system  response  is: 

-  -  »(p)  X(t)  (A-6) 


=  control  torque 
I  =  vehicle  moment  of  Inertia 


P  D(p)  €(t) 


vrtiere 


p  ~  /14- 


d_ 

dt 


T. 


The  initial  conditions  for  the  filter  output  are; 


=  -  (©„  +0  T.) 

0  D  o  A 


€  =  -  0 
o  o 


d^  c 
dt" 


-  0  for  n  >  1 


t=0 


Laplace  transformation  of  equation  (A-6)  then  gives; 


(A-7) 

(A-8) 

(A-9) 


D(b) 


8^  E(s)  -  S  ^  N(p)  ^  (t)j  (A-10) 


Solving  for  E(s); 
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.H(s)  -  V 

8  s'  8  D(b) 


(A-11) 


Tlien: 


€(t)  =  Gp  +  +  O^t  +  G(t) 


(A-12) 


vAiere 


g(t)  ■  iC'V  ^[”(|) ,  A..(t3 

/  0  0(6) 


(A-13) 


If  (t)  in  a  step  function  of  mfpiitudo  X  ,  then: 


g(t)  "  -  A 


-1 


„3 


(A.14) 


Now  swltch-off  will  occur  v/hen  e(t)  -  -  Gp(l-h).  Subni Ltutln,”  this  into 
equation  (A-12)  Gives: 


g(t)  .  -  (hO„  +  9^  .  9^t) 


(A-15) 


The  control  torque  'on-time"  is  equal  to  the  time,  t^,  wliich  satisfies 

equation  (A-15)  for  the  given  value  of  G^,  plus  the  time  delay,  Tlie 

ST/itch  "on-time"  is  t  +  t.. 

o  A 

The  vehicle  rate  at  torque  removal,  G^,  is  given  by: 


t  +  T_ 
o  R 


®f  =  ®o "  / 


(A-16) 
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For  X  (t)  a  atep  functJ.on  of  macnitude 


»f  "  ®0  -  K  (‘o  -  ’n) 


(A-17) 


For  a  given  value  of  can  be  determined  from  equatl'^n  (A-I5)  directly, 

\ 

if  the  expreseion  for  g(t)  is  simply  enougli,  by  iterative  methods  other¬ 
wise,  or  graphically  by  plotting  g(t)  versus  time  and  constricting  the  line 
represented  by  the  right  hand  side  of  equation  (A-15)  -  the  solution  being 
the  intersection  of  the  two  plots.  An  indirect,  but  more  convenient  method, 

however,  is  to  treat  t  as  the  variable.  For  a  given  value  of  t  ,  0  is 

o  0  0 

evalqnted  using  equation  (A-15).  For  this  vflue  of  t^  and  the  correspond- 

*  *  t 

it\^  value  of  0^,  is  determined  from  equation  (A-.l/>).  0^  can  then  be 

plotted  versus  6^,  yielding  the  rate  diagram. 


